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To design a multiferroic material at atomic scale, strong spin-lattice and charge-lattice couplings play crucial 
roles. Our first-principles calculation on (SrCo0 3 ) 1 /(SrTi0 3 ) 1 superlattices, with above coupling properties, 
yields a rich physical phase diagram as a function of epitaxial strain. In particular, a robust ferroelectric 
ferromagnetic insulator of Pc symmetry is stabilized at tensile strain Aa/a 0 = 0.86%-5.53%. The 
polarization can be as large as 36 (iC/cm 2 and magnetic moment can reach 6(i B per unit cell. The 
magnetocrystalline anisotropy energy (0.16 meV/Co in (001) plane, 0.6 meV/Co in (100) plane) is 
comparable with that of TbMn0 3 compound and the magnetoelectric constant a (1.44 X 10~ 3 Gaussian 
unit) is comparable with that of Go 3 B 7 0 13 Br compound. Our study suggests that epitaxially strained 
(SrCo0 3 )i/(SrTi0 3 )i superlattices not only offer an excellent candidate for multiferroic materials, but also 
demonstrate the half-metal and ferromagnetic insulator properties with potential application in spintronic 
devices. 

A ferroelectric ferromagnetic insulator can be a good candidate for multiferroic materials if the cross control 
of polarization or magnetization by magnetic or electric field can be manipulated 1,2 . Thus, a satisfactory 
multiferroic compound not only possesses large polarization and magnetization, but also requires strong 
magnetoelectric coupling 3,4 . As device miniaturization demands integrating electronic and magnetic components 
together, searching for new multiferroic materials constitutes an important research direction 5,6 . 

The perovskite family of compounds offers an excellent candidate for multiferroic materials because of the 
strong coupling among spin, charge, and orbital degrees of freedom 7,8 . In particular, the subtle balance between 
different competing interactions makes the structural tuning very effective so that novel phenomena appear in the 
presence of either compressive or tensile strain 911 . Among others, Lee et al. demonstrated the feasibility of 
converting nonferroelectric and nonferromagnetic EuTi0 3 films into multiferroics under epitaxial tensile strain 
imposed by substrates 12,13 . The critical temperature T c ( = 4.24 K) was very low, it, nevertheless, provided a new 
direction to explore multiferroic materials. The substrate induced strain also helps to make the LuFe0 3 thin films 
a room temperature multiferroic, a weak ferromagnetism was observed at 1 30 K resulting from spin reorientation 
phenomena 14 . Inspired by the experimental advances, the first-principles studies have also offered valuable 
information in this endeavor. The multiferroic property has been predicted by Lee et al. in SrMn0 3 thin films 
under epitaxial strain and T c = 92 K was estimated 15 . An interesting phase transition between a ferromagnetic 
metal and ferroelectric antiferromagnetic (AFM) insulator is also shown in strained SrCo0 3 thin film 16 . There are 
also attempts to transform ferroelectric compounds into multiferroic compounds by doping. For instance, 
ferromagnetism was induced in V-doped ferroelectric La 2 Ti 2 0 7 compound 17 , but the magnetoelectric coupling 
was too weak to function as a multiferroic. Although great processes have been made in the search of multiferroic 
materials, many important issues remain unresolved. The critical temperature is still low which hinders applica- 
tions; compounds with both large polarization and magnetization are rare; even if large polarization and mag- 
netization have been achieved, strong magnetoelectric coupling still needs to be dealt with. 

Since epitaxial strain plays such an important role in tuning material properties, combining crystals with 
similar structures but complementary properties into superlattices can be a good option 18,19 . One expects that the 
resultant materials might combine the characteristics of their constituents. Furthermore, interface mediated 
strain effectively couples the physical properties of their constituents so that the cross controls of constituent 
properties can be realized. In the recent studies on superlattice structures, the improper ferroelectric effect was 
found in PbTi0 3 /SrTi0 3 superlattices 18 and a large magnetoelectric coupling was also observed in nonferro- 
electric NdMn0 3 /SrMn0 3 /LaMn0 3 superlattices 19 . To shed further insight of interface effect on material 
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properties, we have studied, in this Report, the electronic and mag- 
netic properties of (SrCo0 3 ) 1 /(SrTi0 3 ) 1 superlattice (made of 
SrCo0 3 and SrTi0 3 layers of one cubic-cell thick, SCO/STO). As is 
well known, SrCo0 3 is a prototype spin-state transition compound 20 
and possesses strong Jahn-Teller effect while SrTi0 3 n , though 
belongs to paraelectric compound, is near the ferroelectric instability. 
It is hoped that large ferroelectricity and ferromagnetism can arise 
from its constituent compounds while strong magnetoelectric coup- 
ling can be achieved by the interface mediated spin-lattice and 
charge-lattice couplings. 

To explore the phase diagram of SCO/STO superlattices, possible 
crystal structures with different oxygen octahedra tilting (OOT) pat- 
terns are derived from the unstable phonon modes of high symmetry 
P4/mmm structure and their relative stabilities are analyzed. In 
Glazer's notation 21 , the considered metastable crystal structures 
include: P4/mmm (aVc°), P4mm (r 3 ~), Amm2 (r~), P4/mbm 
(a°a°c + ), C2/m (a'aV), Pmna (a^c 0 ), P21/c (a~a-c + ), Pbam 
(P4/mbm with Jahn-Teller distortion), P4bm (P4/mbm + r~), 
Cm (C2/m + r"), Pma2 (Pmna + r~), Pc (P21/c + T~ +r~), 
and Pba2(Pbam + ) symmetries. and refer to the polar 
distortion modes perpendicular and parallel to the (001) plane of 
superlattices (see Supplementary Information). In addition, both 
ferromagnetic (FM) and nearest neighbor antiferromagnetic 
(AFM) states are considered for cobalt moments. Although the 
Co-ion is well known for its varying spin-state depending on crystal 
field, the intermediate spin-state turns out to be the convergent solu- 
tion in the whole epitaxial strain range we studied. In our numerical 
calculations, for a given initial structure with characteristic OOT 
pattern and polar distortion, we search for the convergent metastable 
solution driven by the structural and atomic relaxations. The con- 
vergent solution can be either the one with the same symmetry as the 
initial trial solution or a solution with different symmetry. 

Results 

In Fig. 1, we plot the energies of those possible ground states for the 
strain range e = —5.6% ~ 5.89%. As strain varies from compression 
to tensile, the ground state of the superlattices changes consecutively 
from the antiferromagnetic metal of P4mm symmetry to the fer- 
romagnetic half-metal, insulator, and again half-metal of Pbam 
symmetry, then to the ferromagnetic insulators with almost degen- 
erate P21/c and C2/m symmetries, and is followed by the robust 
multiferroic insulator of Pc symmetry with large polarization and 
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Figure 1 | Energies of SCO/STO superlattices as functions of epitaxial 
strain e for the various possible ground states of different symmetries, e = 

(a — a 0 )/a 0 is defined with respect to the equilibrium lattice constant a 0 = 
3.940A of SrTi0 3 . The vertical lines mark the phase boundaries between 
different symmetries while the dotted vertical line represents the metal- 
insulator boundary within the same structural symmetry. 



magnetization combined with strong magnetoelectric coupling. This 
state is finally replaced by the antiferromagnetic insulator of Pc 
structure. To correlate the electromagnetic properties of possible 
ground states with the specific crystal structures, the lattice constant 
c, OOT, band-gap, and exchange energy are presented in Fig. 2. 
AAFD of Fig. 2b stands for the anti-phase antiferro-distortive vector 
of OOT 22 , its three components are the tilting angles along the 
pseudo-cubic crystal cell ([100], [010], and [001]). Since the AAFD 
vectors for Co0 6 and Ti0 6 octahedra are complementary to each 
other, only AAFD vectors for Co0 6 octahedra are shown. 

Below we discuss the electronic and magnetic properties of various 
possible ground states one by one. For — 5.6% < e £ — 4.3%, the most 
stable state is an AFM metal with P4mm symmetry (see 
Supplementary Information). This state has no OOT and the band 
gaps for both spin-up and spin-down bands are zero. The projected 
electron occupation number and magnetic moment on Co-d orbitals 
are nj°~6.86 and mj°~2.75Li B . The Co ions take an intermediate- 
spin state and form an AFM order. 

As strain e approaches —4.3%, above AFM state is replaced by a 
ferromagnetic state with Pbam symmetry. From the band-gap shown 
in Fig. 2c, this state is initially a half metal with perfect spin polar- 
ization. The half-metal state is further developed into a ferromag- 
netic insulator at e = —3.45%. The band-gaps for both spin-up and 
spin-down bands become nonzero though one decreases while 
another increases with epitaxial strain e. A half-metal state with 
reversed spin polarization could be recovered at e = —0.01% which 
is, however, beyond the border (e = —0.57%) of stability for this 




Epitaxial strain(%) 

Figure 2 | Structural and electronic properties of SCO/STO superlattices 
versus epitaxial strain for the various possible ground states. 

(a) Lattice constant c; (b) AAFD vector. Dashed, dotted, and solid lines 
represent x, y, z components, respectively; (c) Band gaps. Solid and dotted 
lines refer to spin-up and spin-down bands; (d) Exchange energy A£ = 
Eafm ~ E FM . 
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Figure 3 | The polarization of the robust multiferroic state of Pc structure 
as a function of tensile strain. 

state. AAFD vector changes very slowly with epitaxial strain, the 
only nonzero titling angle along z-axis changes from 1.38° to 1.02° 
in the strain range (—4.3%, —0.57%). The electron occupation num- 
ber (n^j 0 —6.90) is almost unchanged while magnetic moment 
(mj°~2.25\i B ) is reduced. These half-metal and ferromagnetic insu- 
lator may find their applications in spintronic devices 23 . 

In the region —0.57% < e £ 0.86% where strain is small, the two 
ferromagnetic insulators with separate P21/c and C2/m symmetries 
appear as the most stable states. These two states are almost degen- 
erate in energy and their crystal structures are all of monoclinic 
crystal class. Similar case was also found previously in EuTi0 3 24 . 
The lattice constants and exchange energies have similar values. 
The electron occupation numbers and magnetic moments are almost 
identical, Wj°~6.90,m^°~2.26u B for both states. The band-gaps for 
spin-up bands are of the same sizes while the band-gaps for spin- 
down differ by a factor of two. This is caused by the different crystal 
structures. For P21/c structure, AAFD vector changes from (2.90°, 
3.10°, 0.02°) to (3.68°, 3.89°, 0°). The x and y components are quite 
close, though not identical due to the orbital ordering similar to that 
in LaMn0 3 23 . For C2/m structure, only x component of AAFD vector 



is nonzero and takes the value 4.17° ~ 5.39° in the corresponding 
strain range. 

At tensile strain e = 0.86%, both energies and first derivatives of 
energies with respect to strain e are identical for P21/c and Pc struc- 
tures. Thus, a second order phase transition takes place from P21/c to 
Pc structure. These two states are quite similar magnetically except 
the polar distortion in Pc structure. The ferromagnetic insulator is 
still kept in Pc structure and the band-gaps evolve continuously 
upwards with e. In the strain range (0.86%, 5.53%) stable for Pc 
structure, the AAFD vector varies from (3.68°, 3.89°, 0.00°) to 
(9.23°, 3.23°, 0.39°). The electron occupation number and magnetic 
moment are almost the same as those of P21/c structure. The total 
magnetic moment is 6/.i B . The most striking feature of this state is the 
large polarization as shown in Fig. 3. The calculated polarization with 
respect to P21/c reference structure lies along [uuv] direction with 
u¥^v. It increases linearly from (0, 0, 0) to ( — 23.9 (iC/cm 2 , 
— 23.9 (J.C/cm 2 , 12.75 uC/cm 2 ) with the inclination angle of about 
27°. This robust multiferroic state is a proper multiferroic state sim- 
ilar to the Pmc2! state found by Lee et al. 15 . It differs from the multi- 
ferroic states studied by other groups where only weak 
ferromagnetism or antiferromagnetism is found to coexist with 
ferroelectric order 14,25,26 . 

Further increasing the tensile strain destabilizes the ferromagnetic 
state and the AFM insulator of Pc structure becomes stable for e > 
5.53%. The x and y components of the AAFD vector are reduced 
significantly while the z component increases considerably. In com- 
parison with the FM state of Pc structure, this state has the similar 
electron occupation number («J°~6.9l) and slightly reduced mag- 
netic moment (m ( f~2.\\\i B s ). This state is also a multiferroic state, 
the calculated polarization is (27.9 uC/cm 2 , 34.0 uC/cm 2 , 32.9 |iC/ 
cm 2 ) at e = 5.89%. 

To shed further insight into the physical properties of the multi- 
ferroic state of Pc structure, the spin and orbital resolved partial 
densities of states (PDOS) are presented in Fig. 4. The epitaxial strain 
is set at e = 2.3% and Fermi energy is set as zero. The Co-d orbitals 
are decomposed into triplet t 2 g{d xy , d yz , d zx ) and doublet 
e g (d x 2_ y 2, d^-r 1 ) sets with respect to Co0 6 octahedron. The two 
nearest neighbors Co-I and Co-II ions are included in order to illus- 
trate charge and orbital ordering phenomena. It is seen that the 
multiferroic state is an insulator with a band-gap of 0.12 eV. 




E(eV) E(eV) 

Figure 4 | The spin and orbital resolved partial densities of states of SCO/STO superlattices at e = 2.3%. The spin-up PDOSs are plotted upwards while 
spin-down PDOS are plotted downwards. Other notations are described in the text. 
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PDOSs are the same for the d, -orbitals of both Co-ions while those 
of d tlg -orbitals are almost the same if x and y-axes are interchanged 
for two Co-ions. To analyze the spin and orbital configurations, the 
projected electron occupation numbers are calculated for Co-ions 
within the muffin-tin radius of 1.25 A. For d t2g -orbitals, all spin-up 
states are almost occupied while only one of the spin-down states is 
occupied for each Co-type. The doubly occupied d tlg -orbital is d zx for 
Co-I and d yz for Co-II. Thus, the Co-ions form an antiferro-orbital 
ordering and take an intermediate spin state. Such orbital ordering is 
made possible by the spontaneous Jahn- Teller distortion similar to 
LaMn0 3 23 (see Supplementary Information). The periodic elonga- 
tion of Co0 6 octahedra along x and y-axes is not only responsible for 
alternately doubly occupied d zx and d^-orbitals, but also responsible 
for the different tilting angles along x and y-axes. We would like to 
emphasize that the Jahn-Teller distortion induced orbital-ordering is 
a typical feature for all the ground states of SCO/STO superlattices. 

So far, we have obtained a ferroelectric ferromagnetic insulator 
with large polarization and magnetization. To verify whether such 
state is a genuine multiferroic state or not, the issues concerning the 
hysteresis loops for ferroelectricity and ferromagnetism as well as 
magnetoelectric coupling have to be addressed. While the hysteresis 
loop for ferroelectricity shows itself automatically by the double well 
structure, the hysteresis loop for the ferromagnetism has to be 
demonstrated through the magnetocrystalline anisotropy energy 
(MAE). By using the algorithms adopted by Xiang et al. 27 , we com- 
pute the MAE of Co-ions of SCO/STO superlattice and the results are 
presented in Fig. 5a-b. The magnetic anisotropy energies are 0.16 
and 0.61 meV at e = 2.3% in the (001) and (100) planes, respectively, 
which have the same order of magnitude as those of Mn-ions in the 
good multiferroic compound TbMn0 3 27 . Because of the oxygen octa- 
hedral tilting, the minima of MAE in (001) plane are located at 35° 
with respect to a-axis. The minimum of MAE in (100) plane is 
slightly off the (001) plane and at 93° measuring from c-axis. 
These data are consistent with the AAFD vector shown in Fig. 2(b) 
where y component (3.65°) is very close to the easy axis. 
Magnetoelectric constant can be similarly estimated using the tech- 
nique devised by Iniguez 28 . The dependences of polarization and 
magnetization as functions of polar-mode distortion are shown in 
Fig. 5c-d. The polar- mode is chosen as a linear interpolation between 



the paraelectric P21/c structure and the ferroelectric Pc structure 17 . 
The dielectric polarity p d n = 55.9|<j| (q is electron charge) and the 
magnetic analogue of the dielectric polarity p™ = 27.64 X 10~ 2 u B /A 
can be extracted from the linear part of P and M. The corresponding 
eigenvalue C„ = 15.07 eV/A 2 is obtained from the total energy 
shown in Fig. 1, and the volume of unit cell is 248. 53A 3 . The esti- 
mated value of magnetoelectric constant is a = 1.44 X 10~ 3 in 
Gaussian unit, which is comparable to that of Co 3 B 7 0 13 Br 29 . To have 
a reliable magnetoelectric constant, as usual, a high convergence 
criterion (10~ 8 eV) is imposed throughout the computations. 
Thus, the tensile strained SCO/STO superlattices do show excellent 
multiferroic property with large polarization and magnetization as 
well as strong magnetoelectric coupling. 

Discussion 

The robustness of the multiferroic state is checked against the dif- 
ferent exchange-correlation functionals and different values of 
Coulomb repulsion U. As shown in the phase diagram for LDA + 
U with U e ff = 3.9eV and GGA + U with U eff = 2.9eV (see 
Supplementary Information), the strain values for transition points 
in Fig. 1 may shift somewhat, but the ordering of various ground 
states are never altered. Thus the overall shape of the phase diagram 
is independent of the specific exchange-correlation functionals and 
the specific value of U. 

The physical mechanism underlying the multiferroic property of 
SCO/STO superlattices can be understood from the crystal structure 
and electron orbital occupation information (see Supplementary 
Information). The asymmetrical bond lengths around Co- and Ti- 
ions in the Pc structure demonstrate off-center polar distortion 
within the Co0 6 and Ti0 6 octahedra. Thus, the electric polarization 
is contributed by both Co- and Ti-ions. The magnetism comes from 
the unpaired Co-d electrons. Since Co-ions are responsible for both 
electric polarization and magnetization, strong magnetoelectric 
coupling is, of course, expected. With regard to the origin of ferro- 
magnetism in the Pc structure, to the zeroth order, the spin-up states 
of Co-ions are all occupied while only one of the spin-down states is 
occupied. The electron occupation is 6 electrons. If we switch to hole 
picture, the level scheme for holes is just the reverse of level scheme 
for electrons. In the hole picture, lowest two Co-d e -orbitals and 
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Figure 5 | The magnetocrystalline anisotropy energy, polarization, and magnetization, (a)-(b) MAE of Co-ions in the (001) and (100) planes, 
respectively; (c)-(d) polar-mode distortion dependence of polarization and magnetization. Solid lines refers to P& M„ dotted lines to P x = P r 
AM X = AM y . Other parameters are described in the text. 
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Co-c?t 2j! -(xy)-orbital are singly occupied by spin- up holes. The nearest 
neighbor Co-d t2g -(yz) and d tlg -(zx) orbitals are alternatively occupied 
by a spin-up hole due to antiphase antiferrodistortive tilting of oxy- 
gen octahetra. The scenario here is quite similar to LaMn0 3 case 
except the exchange of electrons for holes. AFM state would be 
favored if the hole occupation is exactly four. Since the exact electron 
occupation number is large than 6, Co0 2 plane is intrinsically elec- 
tron doped which favors the FM state according to the double 
exchange mechanism. 

We have also extended our study to the generalized (SrCo0 3 ) n / 
(SrTi0 3 ) m superlattices (n = l,2,m = 1,2). The numerical analysis 
suggests that extended superlattice structures do not improve the 
material properties as a candidate for multiferroic compound. A 
close proximity between the ferromagnetic SrCo0 3 atomic layer 
and near ferroelectric SrTi0 3 atomic layer is the best configuration 
ensuring the novel multiferroic property. Deviating from this setting 
will either destabilize the ferromagnetism or favor metallic ground 
state. 

In summary, comprehensive first-principles calculations have 
been carried out for the electronic and magnetic structures of 
SCO/STO superlattices and complex phase diagram has been 
obtained as a function of epitaxial strain. Half-metal and ferromag- 
netic insulator were discovered in the Pbam structure under com- 
pressive strain while a robust multiferroic state was found in tensile 
strained Pc structure. This ferroelectric ferromagnetic insulator has 
large polarization and magnetization, and strong magnetoelectric 
coupling arises from the strong piezoelectric and piezomagnetic 
coefficients as well as the peculiar orbital ordered structure. The 
calculated magnetocrystalline anisotropy energy is comparable to 
that of TbMn0 3 while the estimated interface mediated magneto- 
electric coupling is comparable to that of Co 3 B 7 Oi 3 Br. This robust 
multiferroic state offers an interesting alternative worthy of further 
experimental investigations. 

Methods 

Details of numerical computation. In this Report, the electronic structures were 
studied using a plane wave pseudopotential approach to the density-functional 
theory 30 as implemented in the Vienna ab Initio Simulation Package (VASP) 31 . We 
adopt the generalized gradient approximation with Perdew-Burke-Ernzerhof 
parametrization 32 and the Dudarev implementation 33 for the on-site Coulomb 
repulsion U and exchange interaction J H . The projector augmented wave 
potentials 34 explicitly include 10 valence electrons for Sr (4s 2 4p 6 5s 2 ), 9 for Co 
(3dMs') ; 12 for Ti (3s 2 3p 6 4s 2 3d 2 ), and 6 for O (2s 2 2p 4 ). The wave functions are 
expanded in a plane-waves basis with an energy cutoff of 600 eV. A 5 X 5 X 5 
Monkhorst-Pack k-point sampling is used for a \[l x \fl x 1 (20 atoms) unit cell. 
Each self-consistent electronic calculation is converged to 10~ 6 eV and the 
tolerance force is set to 0.005 eV/A for the ionic relaxation. For Co 3d electrons, 
U e ff = 1.9eV reproduced the ferromagnetic structure with the correct lattice 
constant (a = 3.842A) and magnetic moment (2.967u B ) for cubic SrCo0 3 . This is 
in agreement with previous theoretical calculations 20 and experiment 
measurements 35,36 . For cubic SrTi0 3 , a band insulator with the right lattice 
constant a — 3.940A is also reproduced. To simulate an epitaxial strain, the in- 
plane lattice constant (a = b) of a Vl x \fl x 1 unit cell is fixed while the out-of- 
plane (c) lattice constant and atomic positions are allowed to relax. The 
ferroelectric polarization is computed using the Berry-phase approach 37 . 
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